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Abstract 

The aim of this paper is to present the new results concerning some 
functionals of Brownian motion with drift and present their applications 
in financial mathematics. We find a probabilistic representation of the 
Laplace transform of special functional of geometric Brownian motion us- 
ing the squared Bessel and radial Ornstein-Uhlenbeck processes. Knowing 
the transition density functions of the above we obtain computable formu- 
las for certain expectations of the concerned functional. As an example we 
find the moments of processes representing an asset price in the lognormal 
volatility ans Stein models. We also present links among the geometric 
Brownian motion, the Markov processes studied by Matsumoto and Yor 
and the hyperbolic Bessel processes. 
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1 Introduction 

The aim of this paper is to present the new results concerning some functionals 
of Brownian motion with drift and their applications to financial mathematics. 
The laws of many different functionals of Brownian motion have been stud- 
ied in recent years (see, among others, [5], [3], [7], [T^j, [5D] ), but some of the 
obtained results can not be effective used in application. The distribution of 
Jo du, where = B t + [d with a standard Brownian motion B, is an 
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example of such situation. This distribution can be characterized by Hartman- 
Watson distribution, but the oscillating nature of the last causes the difficulties 
in numerical calculations connected with this functional (see [2]). We study the 
laws of special functionals of geometric Brownian motion, and find results con- 
venient for numerical applications. We investigate the functionals of geometric 
Brownian motion := exp (^B t + fitj for /i 6 R. In particular we study prop- 

y(-l/2) / v 

erties of the functionals T t = d — , , /■>-, — and 1 + BA; for B > 0, where 

:= J*(Yu^) 2 du. We deliver the probabilistic representation for Laplace 
transform of F. In our probabilistic representation of Laplace transform of L we 
use the squared Bessel and radial Ornstein-Uhlenbeck processes. Knowing the 
transition density functions of these processes we obtain computable formulas 
for certain expectations of the concerned functionals. One of the many advan- 
tages of the new result is the fact that they can be effectively used in numerical 
computations. As an example we compute the moments EX", for a > 0, of the 
processes X t representing an asset price in an important stochastic volatility 
model - in the lognormal volatility model. The necessity of computing moments 
results from the problems of pricing derivatives (for instance, the broad class of 
interest rate derivatives necessity a "convexity correction" to the forward rate 
price; for details see e.g. [5]) as well as from the need of approximations of char- 
acteristic functions of random variables with very complicated distributions. 

We now give a detailed plan of this paper. In subsection 2.1 we present 
a method of calculating the moments Er^ for k £ Z (Proposition 12.21 Re- 
marks 12.31 Corollary I2.10p and investigate the connection of functional T with 
a hyperbolic Bessel process (Theorem 12.40 . The general connections between 
hyperbolic Bessel processes and functionals of geometric Brownian motion are 
presented in [T2]. In subsection 2.2 we investigate the different properties of 
1 + . We find two different probabilistic representations of the Laplace 

transform of (1 + BA^)' 1 (Theorems [23] and HT5J}, the form of Eln(l + j3A^) 
and E(l + BA^)^ 1 (Theorem 12.80 . It turns out that for an arbitrary strictly 
positive random variable we can find a representation of the Laplace'a transform 
of (s + for s > 0, in terms of a squared Bessel process (Theorem 12. lip . 
Moreover, we find some interesting connections between E((l + BA^)" 1 ) and 
the conditional expectation of functionals of geometric Brownian motion with 
opposite drift. Notice that we establish all results for a fixed t. Section 3 is an 
illustration of using the previous results in mathematical finance. We assume 
that the asset price process X satisfies dX t = Y t X t dWt with Y being a (GBM) 
(this model is called the lognormal stochastic volatility model or the Hull- White 
model, see [TU]) and Y being an Ornstein-Uhlenbeck process (OU) (the Stein 
model, see [24 J. The distribution of the asset price for the lognormal stochastic 
volatility model is known but degree of complication and numerical obstacles en- 
courage to look for simpler approximations. Jourdain |13| has given conditions 
on existence of the moments, if Y is a (GBM), but not mentioned about how 
to compute it. In this work we find that the moment is equal to the Laplace'a 
transform of the process L (Theorem 13. 1|) . We can also express moments of 
order a > 1 in terms of the hyperbolic Bessel process (Theorem I3.4[) . For the 
model with random time T\ being an exponential random variable independent 
of Brownian motion driving the diffusion Y we find the closed formula of KX^x 
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(Theorem 13. 6p . In Proposition 13.71 we give a closed formula for moments in the 
Stein model. 

Summing up, we present forms of some interesting functionals of Brownian 
motion. Moreover, we find interesting links among a GBM, Markov processes 
arisen during generalization of the so called Pitmann's 2M—X theorem (see |17) . 
|18| . |19) . [201 ) and a hyperbolic Bessel process (see [H]). Finally, we compute 
the moments of the asset price process in the lognormal stochastic volatility and 
Stein models. 



2 Properties of some functionals of geometric Brow- 
nian motion 

Let (ft, J 7 , P) be a complete probability space with filtration F = (J 7 t)tg[o,oo) 
satisfying the usual conditions. Let the process Y be of the form 

Y t = exp(B t -±), (1) 

where B is a Brownian motion. Functionals of Y play a crucial role in many 
problems of modern stochastic analysis. The studies of the properties of integral 
/ Y^du are motivated by the problem of pricing Asian options, (see [15], [2])- 
The process Y t ~ 2 J* Y^du has been considered by Matsumoto and Yor in several 
works concerning laws of Brownian motion functionals. Along with Yf~ 1 J* Y^du 
it plays a central role in a generalization of the Pitmann's 2M — X theorem (for 
details see for instance [T7j, [T5j, [20] )■ Here, we investigate, among others, the 
properties of functional V defined, for /3 > 0, by 

(2) 



1 + PSo Y s ds 



It turns out that this process plays a crucial role in the problem of computing 
the moments of the asset price in the lognormal stochastic volatility model (see 
Section 3). We also find some new properties of the exponential functional 



A W := / (Y^fdu, (3) 
Jo 

where, for 

Y t W :=exp(B t + ^). (4) 

Therefore, Y defined by (JT]) is by definition equal to y^" 1 / 2 ), so T is a func- 
tional of F( -1 / 2 ). We also consider the random variable (which is often called 
a perpetuity in the mathematical finance literature) : 

poo 

4tf := / [Y^fdu. (5) 
We start from investigation of V. 



o 
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2.1 Some properties of T 

Proposition 2.1. // T is given by then Tq = 1 and 

dT t = TtdBt - pT 2 t dt. (6) 

Proof. It follows easily from the Ito lemma. □ 

Proposition 2.2. Let T be given by © andp k {t) = f*ET*du, k eZ,te [0,T]. 
Then the sequence of junctionals (p&) satisfies the following recurrences: 

P'k(t) = 1 + ^ ~ 1] Pk(t) - Pk Pk+1 (t), (7) 

and ^ 

pi(t) = ~E(]n(l + 0jf y u du)). (8) 

Proof. By Proposition 12.11 and the Ito lemma we have 

T k t =l + kf T h u dB u - kp f T k u +1 du + M^Llll T r*d«. (9) 
Jo Jo ^ Jo 



The local martingale J" T k dB u is a true martingale as 

Tl k du < E / K 2fe dw < oo. 



Taking expectation of both sides of Q we obtain . Further 

19-- 



(3dt 



In 1 + /3 / y„du 



o 



as ln(l + /3 / y„du) < ln(l + (3 J Q y u du) and Eln(l + /3 /„ < oo, which 

implies @. □ 

Remark 2.3. Since, by Q, 

Er , = fc(fc_l) pfe(0 _ /3fcpfe+l(t)j (10) 

Proposition 12.21 allows to compute EI^" for k € Z. Taking fe = —1, we easily 
obtain from ([7]) that 

j/_ 1 {t) = l+pt + p- 1 (t), P-x(0) = 0. 

This solution is given by the formula P-xif) = (j3 — l)e* + /3£ + 1 + /3. Notice 
that, having p_i we get recursively from (J7J the functions p_2,p-3, .... Using the 
function pi we can establish p2,P3, ... So, using (|10j) . we can find all moments 
ET^ for fc e Z, provided we know the form of pi. Therefore, to finish this 
computation we need to find the closed form of the function pi . The function 
Pi is given by ([8]), so we have to find E^ln(l + (3 J Q * Y u du)j . The form of pi is 
presented in Corollary 12. 101 
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Now we investigate the connection of T with hyperbolic Bessel processes. 
Let us recall that a diffusion R with the generator given by 

^=-^ + (.+ -)coth(,)-, (11) 

for a E K, is called a hyperbolic Bessel (HB) process with the parameter a (see 
E2 or [3]). Therefore R satisfies 

dR t = dB t + (a + -) coth(i? t )dt. (12) 

We express the Laplace'a transform of functional T in terms of the Laplace'a 
transform of cosh of R. 

Theorem 2.4. Let R be a hyperbolic Bessel process with the parameter a = — 1 
and r be given by (J2J)- For A > we have: 

Ee- XFt =Ee- 0( - cosh{Rt ' ) ~ 1 \ (13) 

where the initial value of the process R satisfies cosh(i?o) = g + !■ 

Proof. Let 6 t = (3T t . Then, by © 

dd t = 9 t dZ t - 2 t dt (14) 

and 6*o = /?. Moreover, for x > 

de"* ' = -e-^xlWt - a;(9 2 dt) + ie-^x 2 ^. 

Taking p(i, x) := Ee~ x6t we get from the last expression that p satisfies the PDE 

with p(0, x) = e _a: ^. Therefore, the Laplace'a transform of t for A > is a 
solution of (fT5|) . Consider a stochastic differential equation (SDE) 



dH t = ^H? + 2H t dB t , H Q = ^>Q. (16) 



Since for any < y < x 



Vx 2 + 2x- vV + 2y < y/(x - y) 2 + 2(x - y) (17) 

there exists a weak solution to SDE ([T6")) and the trajectory uniqueness holds for 
(fig]) (see [H Theorem 5.5.4] and [H Theorem 5.40.1]). Thus, by the Feynman- 
Kac theorem (after changing terminal condition to the initial one in the Cauchy 
problem (fT5"]) ) we obtain that the function u(t,x) := E x e~^ Ht is the unique 
bounded solution of (|15|) withp(0,x) = e~ x ^ (see |14l Theorem 5.7.6]). Let us 
define the diffusion St := i?t + 1. It is easy to check that 

dSt = s[sJ^ldB t , H = j 3 +l. (18) 
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By the same arguments as before there exists a weak solution to (fT5|) and the 
trajectory uniqueness holds for (fT5)) . Now, we observe that the diffusion J7 t = 
cosh(i? t ), where R is the hyperbolic Bessel process with the parameter —1, and 
such that cosh(i?o) = a + 1 is the solution of (Tig)) . So, the processes S and U 
have the same law. Thus, 

Ee- Art = Ee-? e * = p(t,X//3) = Ee~ fiBt = e p Ee~ ps * = e^E e -^ cosh(Rt K 

This ends the proof. □ 

The connections between hyperbolic Bessel processes and functionals of ge- 
ometric Brownian motion are presented in |12| . 



2.2 Some properties of (1 + fiA^) 

We start form the computation of the Laplace'a transform of (1 + jUA^' 
It is worth to remark that we compute it for a fixed time t. We can find in 
literature (see for instance [T^]) that the problem of computing of expectations 
for functionals of geometric Brownian motion for a fixed time is in general much 
difficult than with stochastic one (see also Subsection 13. 2. 21) . 

Let us recall that a squared ^-dimensional radial Ornstein-Uhlenbeck process 
with the parameter —A for 6 > 0, A € M, is the solution of the SDE 

X t = x+ [ {6 -2\X s )ds + 2 [ y/X~ s dW s , (19) 
Jo Jo 

where W is a standard Brownian motion. For detailed studies of these processes 

see [5] and [3]. If A = 0, then the strong solution of (|19p is a squared 6- 

dimensional Bessel process (see [H]). The number 6/2 — 1 is called the index 

of the process. In the sequel we will use the notation X x for the process X 

starting from x, i.e. Xq = x. 

In the next two theorems we find a probabilistic representation of the Laplace 

transform of (1 + fiA^y 1 . 

Theorem 2.5. Assume (3 E (0, 1], /jgl and t > 0. Then, for any A > 

( ^)) = E M0 X (- ln(v^))), (20) 

v 1 + BAT } ] 



i exp 



where Q x {t) is a squared 0- dimensional radial Ornstein-Uhlenbeck process with 
the parameter —1 such that 6(0) = A and, for x > 0, 

&(a:)=Vt(l,aO, (21) 

and for x > 0, s > 

Ms,x)=EG t (R x (s/2)). (22) 
Here R x is a squared Bessel process with the index —1 starting from x, and 

G t (x) = e^ 2 / 2 Ecxp (jjiB t + ^(s t 2 - ^(B t ))), (23) 
B is a standard Brownian motion, and 



f x{y) = In (xe y + cosh(y) + \J x 2 e~ 2y + sinh 2 (y) + 2xe~ y cosh(y) 
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Proof. Let B be a standard Brownian motion under P. Let us define the function 



:=Ee^- — _pj (24) 

for x > 0, s G [0,1]. Observe that q(s,x) < 1. It is not difficult to check, 
using the Lebesgue theorem, that q belongs to the class C 1,2 ([0, oo) x [0,oo)). 
Moreover, it is easy to see that q satisfies the partial differential equation 

- a ^l =x (^l + ^l) ) (25) 
ds V dx dx 2 / ' 

and q(s,0) = 1. Define 

0t(x) := <?(l,a;) = Eexp ( ~~TJTj") and P^^) := q(e~ s , a;). 

1 + A t 

Then p(s, x) satisfies the partial differential equation 

(26) 



dp / <9p dp \ 
9s V 9a; dx 2 J ' 



s > 0, x > and p(0, x) = 4>t(x). Consider a diffusion {/ with a generator of the 
form 

d 2 d 

A u = x -n + x ^-- ( 27 

dx z ax 

This diffusion satisfies the SDE 

dU t = V2^/u~ t dW t + U t dt, (28) 

where W is a standard Brownian motion. Since p(s,x) < landpG C 1,2 ([0,oo)x 
[0,oo)), using the Feynman-Kac theorem (after changing terminal condition 
to the initial one in the Cauchy problem (|2l)|)) we obtain that p admits the 
stochastic representation p(s, x ) = K(f>t(Ug) (see |141 Theorem 5.7.6]). 

Observe that 9(s) := U2 S is a 0-dimensional radial Ornstein-Uhlenbeck pro- 
cess with the parameter —1, as ([28]) takes, by the scaling property of Brownian 
motion, the form 

d0(t) = 2y/6(t)dB t + 26(t)dt. (29) 

Thus 

q{s,x) = p(-lns,x) = W4 t (6 x (-(l/2)lns)) = E<f> t {6 x {-]n^)), (30) 
9(0) = x. Hence, taking x = X, s = f3 we see that (|30p gives ([2U|: 

Eexp( ^—-) =q(s,X) =EM0 X (-HVP)))- 

v 1 + l3Af> 1 

To finish the proof we have to find the form of <j> t . To do this we define the new 
functions: 

■4>t(s,x) := Eexp ( ^777). (31) 

v s + A\ M ' 

G t (x) :=Ecxp(-^ T ), (32) 
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for s > 0, x > 0. So ijj t (l,\) = 0t(A) and Gt(x) = ipt(0,x). Observe that ip t 
satisfies the partial differential equation 

ds - cb 2 • 1 ' 

Consider a diffusion X with a generator of the form 

Ax = X dx^- 

Using again the Feynman-Kac theorem we deduce that ip admits the stochastic 
representation ip t (s,x) = EGt(Xf). Now observe that R(s) := X% s satisfies 

dR(t) = 2y/RjrjdB tl (34) 



so R is a squared Bessel process with the index — f . Therefore, we obtain 
It remains to compute the form of the function Gt- Define a new probability 
measure Q by 

§U --»(-<*-£); < 35 > 

Since B is a standard Brownian motion under P. then B t = B t + fJ,t is a standard 
Brownian motion under Q, by the Girsanov theorem. For Af^ := e 2 ® u du we 
have 4 M) = ^ 0) ; so 

G t (x) =E e ~^ =E Q f e "^ T e^ + ^ tN ) =e^ 2 */ 2 E Q exp f ^ 



(0) 



A 

Now we use the Matsumoto-Yor result |I9l Thm. 5.6], which states that 
where 

fx{y) = argcosh(cce _? ' + cosh(y)) 



= In ^x,e y + cosh(y) + y x 2 e~ 2y + sinh (y) + 2xe~ y cosh(y) 
In result we obtain 

G t (x) = e-^ 2 / 2 Ecxp (jjiB t + ^(b 2 - ^(fl t ))), 
where B is a standard Brownian motion under P. This finishes the proof. □ 

Our next theorem provides another probabilistic representation of Laplace 
transform of (f + pA[ now for /3 > 0. 

Theorem 2.6. Fix (3 > 0, fi € R and t > 0. Then, for any A > 0, 

Ecxp ( ^) - EG t (i? A ^(l/(2/3)), (36) 

where R x ^ is a squared Bessel process with the index A//3 starting from x and 
Gt is defined by 
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Proof. Formula follows from the proof of Theorem l2.51 since for s > 0, x > 
we have, by and ([22]). 

Eexp( =^(«,x) =EG t (i2 s (»/2)) (37) 

and taking x = X/(3 and s = 1//3 we obtain ([M]). □ 
Corollary 2.7. for /3 G (0, 1] we have 

EM0\- ln(v^))) = EG t (i? A ^(l/(2/3)), (38) 
where (9 A is defined in Theorem \2.5\ 

Proof. (|3S) follows from © and J2DI). □ 
Using this result we can obtain the expectations of ln(l + (3A[^) and (1 + 

Theorem 2.8. Fix (3 > 0, fx e K and t > 0. Then 



Eln(l + /3A w ) = / -^^(l-e-^)dy, (39) 

Jo y 

E ( ^l) = 1 - P r Gt(y)e- yfi dy, (40) 

where Gt is given by (|23p . 

Proof. Let f(f3) = Eln(l + / 8A t (,i) ) for /3 > 0. Since 

E| ln(l + l3A { / } )\ < 1 + (3EA { t ^ < oo, 
the function / is well defined, continuous and /(0) = 0. Moreover, for j3 > 0, 

A^ \ 1/, _/ 1 



/'(/3) = E — — - ) = — ( 1 — E( )). (41) 

By definition of ip (see (|3"Tj) ) and ([2"2"]) we know that for s > 0, a; > we have 

Eexp( ^) = Tp t (s,x)=®G t (R x (s/2)), (42) 

s + A t 

where R x is a Bessel process with the index —1 starting from x. Since, by 
definition (see ([3^]) ). G*(0) = 1 and the transition density functions for the 
process R x are known (see [2T] Chapter IX, Corollary 1.4]) we can write 

EG t (R x (s/2)) = e- x ' s + [°° G t (y)-, f^e-^/ s I 1 (2^Ey/s)dy, (43) 

Jo s V V 
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where I\ is the modified Bessel function. Let us recall that (Ii(x)/x)' = l2(x)/x 
(see [H Appendix 2]). Hence and by (g3]) and (gj) we obtain, for x € [0, 1], 

-E( 1 ^ exp f X Vk- 



x/s+ i(.r Gt(y) j^ ( ^ )/s/i(2 ^ /s) ^ 

/•OO 

(1 - J)( J G t (y)e-y/ s (2^/ S y 1 I 1 (2^/s)dy 



1 -«/, 2e~ x/ 
—e x/s H =- 



s i 

+ — a"! / ^ G t {y)e-y/ s {2^/ S )- 1 I 2 {2^/s)^dy 
s K Jo s 

Let x tend to 0. Since we can pass with the limit under the integrals, using the 
asymptotic behavior of Bessel functions, i.e. Ii(x)/x ±2 1/2 and li{x)jx ±2 x/8 
in neighborhood of (see [JJ Appendix 2]) we obtain 

/ 1 \ 1 1 f°° 



or equivalently 



s 

Putting s = 1/(3 in (g5J) yields (gDJ). From (jUJ and (gDJ) we conclude 

/'(£) = / G t (y)e-y^dy. 



This finishes proof of the theorem, since /(0) = 0. □ 

Remark 2.9. Formula (gDJl gives the closed expression of E((l + PAp*)- 1 ) 
for P > 0. The density of A^ is known in literature, but due to complicated 
nature of Hartman- Watson distribution, it can hardly be used for numerical 
computations (see for instance [19] and [2]). Since the simple form of function 
Gt is given explicitly, the formulae (|59")1 and (|4H|) allows to obtain numerically 
Eln(l + pA^) and E((l + pA^)" 1 ). 

Theorem 12.81 allows to find the first function pi(-) for the recurrence estab- 
lished in Proposition 12.21 

Corollary 2.10. Let pi be given by ©. Then 



Pl (t)4f^M(l-e^ 

P Jo y 



)dy, (46) 

P Jo V 

where G is defined by (|23p . 

Proof. Since Z u = (1/2)B± U is a standard Brownian motion, we infer that 

it pt 



□ 



Pl {At) = Eln(l + /3 / Y u du) =Eln(l + 4/3 / e^ 4 "^"^) 
Jo Jo 

= Eln(l + 4/3 / e 2(z " - "W) = Eln(l + 4^4 _1) )- 
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Let us now observe that we can deduce more general fact from the proof of 
Theorems l2.5l and l2.8l It turns out that for an arbitrary strictly positive random 
variable we can find a representation of the Laplace'a transform of (s + £) _1 for 
s > 0, in terms of a squared Bessel process. It gives also a second simple proof 
of general version of (J22J) • 

Theorem 2.11. Let £ be a strictly positive random variable. Then for any 
x > 0, s > 

A =EG(iT(s/2)), (47) 



-exp 



where R x is a squared Bessel process with the index —1 starting from x and 



G{x) = Eexp 



Moreover for /3 > 



:in(i + /?o 



G(y) 



(i 



Idy. 



(48) 



(49) 



Proof. Let's take a copy of R x independent of £. Then (|48|) implies 

R x {s/2Y 



EG(R x (s/2)) =Eexp 

= Eexp I - 



l+e^s 



EE exp -C R x {s/2] 



-exp 



where we used the form of Laplace transform of squared Bessel process (see |211 
Chapter XI, page 441])). The proof of (j4U|) goes in the same way as in Theorem 
HE □ 

Now, we use formula (|4"0)) and the results of Matsumoto and Yor to ob- 
tain some interesting connections between E((l + fiA^) -1 ) and the conditional 
expectation of functionals of geometric Brownian motion with opposite drift. 

Proposition 2.12. For /i > and j3 > we have 

1 



E 



1 + 2/3,4 



Cm) 



l-2/m(4-' i) \A£ti = 1/(2/3) 



(50) 



Proof. By the result of Matsumoto and Yor [THl Thm. 2.2] the process 



0} on the set {A, 



{B 



> 



1/(2/3)} has the same distribution as the process 



log(l + 2/3 A ( f'),t > 0} for fi > 0. From that we obtain 



E 



(4^1^ = 1/(2/3) 



E 



o (l + 2/3A i f ) ) 2 



ds = — (l-E( —r^ 



Proposition 2.13. Lei /3 > and /z £ 

E(Ai~^\A^ = 1/(2/3) 



Then 
1 

" 2 



GtWe-vPdy. 



□ 



(51) 
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Proof. It follows from Proposition 12.121 and gOj) . □ 
Proposition 2.14. For (3 > 0,/i > we have 



2/iB ( ~' j) 

K 'w-J = 1 - 2 ^(^""VL-^ - V(2/3)). (52) 

Proof. Fix /i > 0. Define the new probability measure Q by 

d® = 2MBt -a M » t (53) 

The process Vt = B t + 1[it is a standard Brownian motion under Q, so 

/ 1 x / 2 M (v t -Mt) >. , ea/iB'-"' 

E( ^ = Eq ( 1 = E( ? — r), 

Vl + 2/M^' H + 2/3j 'e 2 / "-'«W M + 

Now the thesis follows from Theorem 12. 121 □ 

Proposition 2.15. For t > we have 

Pi(t) = t-4pJ E(4- 1) |A^ 1 )=l/(4^))ds. (54) 

Proo/. We have = ET t , pi(0) = 0, and 

r « = =7317 > ( 55 ) 

1 + 4/3A] ' 

where B t = B± t /'2 is a standard Brownian motion, and ^ is defined by 
with i? instead of B. Since 

E^V^ 15 = V(4/3)) = E^Vo^ = 1/(4/3)), 

the thesis follows from (|5"2")) with /u = 1. □ 

Remark 2.16. Notice, that we establish all the results for fixed t. In several 
papers (for instance |T3] , HH]) the integral functionals of a geometric Brownian 
motion with random time given by random variable independent of Brownian 
motion and with exponential distribution were investigated. In particular, 



Eln (l + /3 j Y^duj = Eln (l + /3^p), (56) 

because f^ x Y 2 du — d where Ci.a is a random variable with beta distri- 

bution with the parameters 1 and a = V 2X+1 ^ 4 1 / 2 ^ ^ b jg a ranc i om variable 

with gamma distribution with the parameter b = V 2X+1 ^ 4+1 / 2 ^ ^ ^ an j ^ are 
independent (see |19|). Later, we also explore idea of using random time. In 
subsection [322] we show how to compute the moments in a lognormal stochastic 
volatility model with random time being exponentially distributed and indepen- 
dent of Brownian motion driving the model. 
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3 Moments of the asset price in the lognormal 
stochastic volatility and Stein models 



3.1 Model of market 

We consider a market defined on a complete probability space (0, J 7 , P) with 
filtration F = (J r t ) t6 [o,T]j T < oo, satisfying the usual conditions and J- = Tt- 
Without loss of generality we assume the savings account to be constant and 
identically equal to one. Moreover, we assume that the price X t of the under- 
lying asset at time t has a stochastic volatility Y t being a geometric Brownian 
motion or an Ornstein-Uhlenbeck process, so the dynamics of the proces X is 
given by 

dXt = Y t X t dWt„ (57) 

where Xq = 1. In case of Y being a GBM the dynamics of the vector {X, Y) is 
given by (|5?j) and 

dY t =Y t dZ u Y = l, (58) 

and in case of Y being an OU the dynamics of the vector (X, Y) is given by 
([STjl and 

dY t = -XY t dt + dZ t , Y = l (59) 

for A > 0. The processes W, Z are correlated Brownian motions, d(W, Z) t = pdt 
with p G [—1,1]. In the both cases the process X has the form 

X t = e fo Y ^-fo Y ^/^ (60) 

and this is a unique strong solution of SDE (|57p on [0, T]. The existence and 
uniqueness follow directly from the assumptions on Y and the well known prop- 
erties of stochastic exponent (see, e.g., Revuz and Yor [H]). Since the process 
X is a local martingale, there is no arbitrage on the market so defined. 
Notice that we can represent W as 

W t = pZ t + Vl - P 2 V t , (61) 

where (V, Z) is a standard two-dimensional Wiener process. Using and (|6"Tj) 
we can expressed the moment of order a of X as 

EX? = Ee Qp ^ Y u dz u+a ^T^^ Y u dV u -fft Y ' du (62) 

= Ee ap ti Y u dz u + " 2 ^-.f ) - a So y*du_ 

3.2 Moments of the asset price in the lognormal stochastic 
volatility model 

3.2.1 Moments of order a > 

In this subsection we will calculate moments of order a > in the lognormal 
stochastic volatility model, so for Y of the form ([55)) . Jourdain [T5] gave a 
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sufficient condition on a > 1 for existing of moments of X. Namely, Jourdain 
proved that for a 6 (1,(1 — p 2 ) -1 ) and p ^ the moments EX a exist, but he 
didn't find the value of these moments. We calculate the value of moments for 
a > 0, a(l — p 2 ) < 1 and p € [—1, 1]. Sin [231 established that the process X 
is a true martingale if and only if p < 0. First, we prove that the moment of 
order a of the strong solution of (|5T|) is equal to the Laplace'a transform of the 
process T. 

Theorem 3.1. Let t € [0,T], a > 0, a(l - p 2 ) < 1 and T 6e #roen by ©. J/ 
X is given by (|57p . £/ien 

EX t Q = e - {/3+ " Q) E exp ((/? + pa)r t ) , (63) 

w/iere 



/8 = ^/a-a 2 (l-p 2 ). (64) 
Proof. Define a measure Q by 

^ =e -p j? tfy:*", (65) 

C^P .7-V 

where /3 is given by (|64p . The measure Q is a probability measure since, by 

USE" 



Using (|62p and the definition of Q we infer 

EX" = Ee a »ti Y u dz u +^p^ Po Y-du = Eqe ( P a+p)(Y t -i) . (66) 

By the Girsanov theorem, B t = Z t + f j3Y s ds is a standard Brownian motion 
under Q and 

dZ t = dB t - f3Y t dt, Z a = 0. (67) 

We know, by the result of Alili, Matsumoto and Shiraishi [TJ Lemma 3.1], that 
the unique strong solution of (|67]l is given by 
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1 + P Slu s ds^ 



where 
Therefore 



1 + P J* U s ds 



The law of the process Y under Q is equal to the law of the process T under 
P, since the law of the process Y under P is equal to the law of the process U 
under Q. Hence and by we obtain (|S3"]l . This ends the proof. □ 



15 



Remark 3.2. a) From Theorem [XT] we immediately see that all moments exist 
provided p 2 = 1. 

b) The condition a(l — p 2 ) < 1 is not a necessary condition for existence of 
moments since in case of p = the process X is a martingale, so EX t exists. 
Although, for p = 0, EX t a = oo for a > 1 ([13]). 

Remark 3.3. In Theorem 13. II we prove that the computation of EX" for a > 
such that a(l — p 2 ) < 1 is equivalent to the computation of the Laplace'a 
transform of T t at point A = (3+ pa. In turn, the recurrence from Proposition ^. 21 
allows to find EP^ , so we can find an approximation of the Laplace'a transform 

of r in the neighborhood of zero by its moments, namely Ee Ar * w J2iLo ^~7r 
for sufficiently large N. In this way we obtain an approximate value of EX" . 

We can also express moments of order a > 1 in terms of the hyperbolic 
Bessel process with the parameter — 1. 

Theorem 3.4. Assume that a > 1, a(l — p 2 ) < 1. Let X be given by (|5T|) 
with Y given by (JTJ) and R be a hyperbolic Bessel process with the parameter 
— 1 . Then 

EAT t a = e- pa Ee-P cosh ( Rt \ (68) 
where cosh(i?o) = —pa./ '(3, fi = \J a — a 2 (l — p 2 ). 

Proof. We use Theorem 13.11 and Theorem 12.41 with A = — (f3 + pa), and A > 
provided a > 1. □ 



3.2.2 Moments with independent random time 

In this subsection we find the closed formulae for the moments in a lognormal 
stochastic volatility, when the time is an exponential random variable indepen- 
dent of Brownian motion driving the diffusion Y. The idea of considering such 
a time is not new and can be find in many studies of Asian options (see for 
instance (TS], |19|). 

Proposition 3.5. Let T\ be a random variable with exponential distribution 
with the parameter A > 0. Assume that T\ is independent of a standard Brow- 
nian motion B. Let Z t = 2Bt, Y t = e~? +Zt and U t = e Bt ~~ f . Then 



E(ln(l + pj Y u duj} = y - 4/3 2 J A j U 2 ds - K/4) + (l + PL\)- 2 dK, 

(69) 



E( / U 2 ds - K/4) + = r / e~ u u ^ \ l-Ku/2) du. 



o XT 



V2A+1-1 
2 



(70) 



Proof. It is obvious that Y 4t = e ~ 2t+Zlt = e ~ 2t+2Bt = Uf . Using the Taylor 
theorem with integral remainder to the function f(x) = ln(l + /3x) gives 

/•OO 

ln(l + fix) = (3x - S3 2 \ {x- K) + {1 + l3K)~ 2 dK. (71) 



16 



Hence replacing x by Jq Tx Y u du and taking expectation we get 



E(ki(l + pj Y u du)^ = /3e( j Y u duj (72) 

poo <-4T A 

- P 2 \ E( / K u du - 1T) + (1 + PK)~ 2 dK 
Jo Jo 

= 4/3E( / A [/ 2 ds) - 4^ 2 / E( / A E7?ds - tf/4)+(l + 0K)~ 2 dK. 



Let A t = / Q * U 2 ds. The Mansuy and Yor theorem [E2 Thm. 6.1] gives £70]) and 
EA Ta = 1/A. This and (JZ2J) completes the proof. □ 

In the next theorem we establish the explicit formula for moments of Xt 2X ■ 

Theorem 3.6. Let a > 0, a(l — p 2 ) < 1 and T\ be a random variable with ex- 
ponential distribution with the parameter A > 0. Assume that T\ is independent 
of Brownian motions V and Z driving the process X. Then 



2Tx A r(i + ViATT) v 1 

1-1/2/3 <-oo _ 

U x (1/2/3) / e s ^-y(h(y)dy + 02(1/2/3) / e^^i(y)dy 

V JO Jl/2/3 



where {3 = a — a 2 (l — p 2 ), 

= a --( 1 +^+4A)/2$^( 1 + VTT4A)/2, 1 + VTTlA,^- 1 ), 

2 (.t) = x- (1+ ^ TW)/2 *((i + VTTlA)/2, i + VTTIa,^- 1 ), 

and $, ^ denote the confluent hypergeometric functions of the first and second 
kind, respectively 

t / . \ - (a) k z K 



k=0 



iYl-V) Th — 1) , 

tt(a, 7, «) = rn V_ 7J 7, *) + V, , V ^S 1 + a - 7) 2 - 7,0), 

r(l + a — 7) r(a) 

where (a)o = 1 emd 

, , r(a + fc) . . 
r(a) 

fork= 1,2,... 

Proo/. If = ±Z 4t , 5 t = e 5 *-*, then K 4t = S" 2 . So 

dS t = S t (dB t - idt), So = 1, 

and 



1 + /3 J Y s ds 1 + 4(3 f* S 2 du 



(74) 
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Define a new process 

_ 1 l + wfifidu 

From the Ito lemma 

d9 t = -26 t dB t + (40* + l)dt, O 

Now we observe that the diffusion 9 has the generator 

d 2 . d 

+ 4z + l 
dx* dx 

which is identical with the generator of the process 



(75) 



1 

w 



Ae = 2x 2 — + (4x + l) — 1 (76) 



Xt = cxp(2S t + 2t)(^ + J cxp(-2B u - 2u)duj 
since from the Ito lemma 

d Xt = 2xtdB t + (4 X t + 

Hence and from the fact that xo = &o we deduce that processes 9 and x have 
the same distribution. Let us take another Brownian motion BT = \FlBt and 

1 2 

define the process 



exp(V2B; + + jf cxp(-V2i?: - u)d«) . (77) 



Using the fact that 9 t and x* have the same distribution, we obtain 29 1 = d r/ t . 
Moreover, we know that ry t is a Markov process with the resolvent 



□ 



Uxf(x) = + V4 3±4( 2) Ui{x) I* e-^ 2 (y)f(y)dy + Mx) [°° e-^My)f(y)dy). 

1(1 + V4A+1) v Jo Jx ' 

(78) 

(for details see [7J Theorem 3.1]), so we conclude by Theorem 13.11 (fTI)) . 
and definition of r\ that 

= e-^Ee X p{^-L} = ie-^^(i), 
with/(a;)=exp{^i}. 

3.3 Moments of the asset price in the Stein model 

In this section we consider the Stein model, i.e. the model described by (|57p 
with Y being an Ornstein-Uhlenbeck process, so Y is given by ([M]) : 

dY t = -\Y t dt + dZ u Y = 1, A > 0. 

For t in neighborhood of zero we find an exact value of EX". Let b be the 
unique solution of equation 

6(1 - e - 2b ) = 2. (79) 
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Proposition 3.7. Let a > and p be such that a(l — p 2 ) < 1, p < X/a and 

7 2 = A 2 -a 2 (l - p 2 ) + a- 2\pa > 0. (80) 

Ift G [0,|), i/ien 

2/3 

where b is given by (f?9")) and 



EXf = e (1+t) ' 3 (cosh( 7 i) + — sinh( 7 i)) 2 , (81) 



j 9 = i(A-pa). 

Proof. Define the new measure 



r/P 



Clearly, Y" t is the Gaussian random variable with the mean e and variance 
i(l - e~ 2At ). Moreover EeT~ F « du < oo for u < t, since ^(1 - e~ 2At ) < 1 by 
assumption on t. Thus, by Jensen inequality, 



t Jo ' t Jo 

^ ft e~ 2A "«A 2 /(2-tiA(l-cxp(-2Ati))) 

du 

t 



o 



1 - 3f (1 - e- 2A ") 



1 f* e uX2 / 2 
< — I — = du < oo. 



l-^(l-e- 2A ") 



So, Q is a probability measure, by the Novikov criterion. Observe, by the 
Girsanov theorem, that the process Y is a Brownian motion under Q starting 
from 1. Formula ([59")) implies 

J Y u dZ u = i(F t 2 — (t + l))+X J Y 2 du, 
so by we obtain 

/, M \afl-A aa-X v-2 , 2Aap-A 2 +a 2 (l-p 2 )-cv r t v 2 j 

EX t a = e- {t+1) -^-E Q e-^- Y t + * 5° Y « du 

= eS 1+t)l3 Ee-' 3B t B ^ du 

= e™> ( cosh (7 t) + M sinh(7 t)) "* exp (t ( ^ /2+ 2 f 7t 

V V ; 7 W V V 2 C osh( 7 i) + M S inh( 7 t) 

where in the last equation we use the form of the Laplace'a transform of [B 2 , J ' B 2 du), 
where B is a Brownian motion starting from 1 (see e.g. [H formula 1.9.7 page 
168]). □ 

Remark 3.8. If a € (0, 1), then flSOj! holds. If a > 1, a(l - p 2 ) < 1 - 2Ap and 
p < then flSE) holds. 
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